A supermassive scalar star at the Galactic Center? 
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ABSTRACT 

We explore whether supermassive non-baryonic stars (in particular boson, mini-boson and 
non-topological soliton stars) might be at the center of some galaxies, with special attention to 
the Milky Way. We analyze, from a dynamical point of view, what current observational data 
show, concluding that they are compatible with a single supermassive object without requiring 
it to be a black hole. Particularly, we show that scalar stars fit very well into these dynamical 
requirements. The parameters of different models of scalar stars necessary to reproduce the 
inferred central mass are derived, and the possible existence of boson particles with the adequate 
range of masses is commented. Accretion to boson stars is also analyzed, and a comparison with 
another non-baryonic candidate, a massive neutrino ball, which is also claimed as an alternative 
to the central black hole, is given. Both models are capable to explain the nature of the object 
in Sgr A* without invoking the presence of a singularity. One difficult issue is why the accreted 
materials will not finally produce, in a sufficiently long time, a black hole. We provide an 
answer based on stellar disruption in the case of boson stars, and comment several suggestions 
for its possible solution in neutrino ball scenarios. Finally, we discuss the prospects for the 
observational detection of these supermassive scalar objects, using the new generation of X-ray 
and radio interferometry satellites. 
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During the last years, the possible existence of a 
single large mass in the Galactic Center has been 
favored as the upper bound on its size tightens, 
and stability criteria rule out complex clusters. Al- 
though it is commonly believed that this central 
mass is a supermassive black hole, it is not yet 
established, as we discuss below, on a firm obser- 
vational basis. 

The aim of this paper is to present an alterna- 



tive model for the supermassive dark object in the 
center of our Galaxy, formed by self-gravitating 
non-baryonic matter composed by bosons. This 
kind of objects, so-called boson stars, are well 
known to physicists, but up to now, observational 
astrophysical consequences have hardly been ex- 
plored. The main characteristics of this model arc: 

1. It is highly relativistic, with a size com- 
parable to (but slightly larger than) the 
Schwarzschild radius of a black hole of equal 
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2. It has neither an event horizon nor a sin- 
gularity, and after a physical radius is 
reached, the mass distribution exponentially 
decreases. 

3. The particles that form the object interact 
between each other only gravitationally, in 
such a way that there is no solid surface to 
which falling particles can collide. 

It is the purpose of this work to show that 
these features are able to produce a Galactic Cen- 
ter model which can be confronted with all known 
observational constraints, and to point the way in 
which such a center could be differentiated from a 
usual supermassive black hole. 

1.1. Why scalar fields? 

Interesting models for dark matter use weakly 
interacting bosons, and primordial nucleosynthesis 
show that most of the mass in the universe should 
be non-baryonic if SI ~ 1. Most models of inflation 
make use of scalar fields. Scalar-tensor gravitation 
is the most interesting alternative to general rel- 
ativity. Recent results from supernovae, which in 
principle were thought of to favor a cosmological 
constant (Perlmutter et al. 1999), can as well be 
supported by a variety of models, some of them 
with scalar fields too (Celerier 1999). Particle 
physicist expect to detect the scalar Higgs particle 
in the next generation of accelerators. Scalar dila- 
tons appear in low energy unified theories, where 
the tensor field g^ v of gravity is accompanied by 
one or several scalar fields, and in string effective 
super-gravity. The axion is a scalar with a long 
history as a dark matter candidate, and Goldstonc 
bosons have also already inferred masses. Sym- 
metry arguments, which once led to the concept 
of neutron stars, may force to ask whether there 
could be stellar structures made up of bosons in- 
stead of fermions. 

In recent works, Schunck and Liddle (1997), 
Schunck and Torres (2000), and Capozziello, Lam- 
biase and Torres (2000) analyzed some of the ob- 
servational properties of boson stars, and found 
them notoriously similar to black holes. In par- 
ticular, the gravitational redshift of the radiation 
emitted within a boson star potential, and the ro- 
tational curves of accreted particles, were studied 
to assess possible boson star detection. The in- 
terest in observational properties of boson stars 



also led to investigate them as a possible lens in 
a gravitational lensing configuration (Dabrowski 
and Schunck 2000). Recent studies are putting 
forward the gravitational lensing phenomenon in 
strong field regimes (Virbhadra et al. 1998; Virb- 
hadra and Ellis 2000; Torres et al. 1998a,b). This 
would be the case for boson stars, which are gen- 
uine relativistic objects. 

This work is organized as follows. Sec. II is a 
brief summary of the main observational results 
concerning Sgr A*, and the main hypotheses in 
order to explain it are outlined there. Sec. Ill an- 
alyzes what dynamical observational data show, 
and what kind of models can support them. Sec. 
IV gives the basic ingredients to theoretically con- 
struct scalar stars, shows mass and radius esti- 
mates, and study effective potentials and orbits 
of particles. In Sec. V we study the center of the 
Galaxy, show which are the stable scalar star mod- 
els able to fit such a huge mass, and comment on 
the possible existence of boson particles with the 
required features. We also provide there an as- 
sessment of the disruption processes in boson star 
scenarios. Discussion and conclusions are given in 
Sections VI and VII. 

2. The Galactic Center 

2.1. Main observational facts 

The Galactic Center is a very active region to- 
ward the Sagittarius constellation where, at least, 
six very energetic radio sources are present (Sgr 
A, A*, Bl, B2, C, and D). Furthermore, there are 
several supernova remnants, filaments, and very 
reach star clusters. Several observational cam- 
paigns (Gcnzcl et al. 1994) have identified the ex- 
act center with the supermassive compact dark 
object in Sgr A*, an extremely loud radio source. 
The mass and the size of the object has been estab- 
lished to be (2.61 ± 0.76) x 10 6 M Q concentrated 
within a radius of 0.016 pc (about 30 lds)(Ghez 
1998; Genzel et al. 1996). 

More precisely, Ghez et al. (1998) have made 
extremely accurate velocity measurements in the 
central square arcsecs. From this bulk of data, it 
is possible to state that a supermassive compact 
dark object is present at the Galactic Center. It 
is revealed by the motion of stars moving within 
projected distances around 0.01 pc from the radio 
source Sgr A*, at projected velocities in excess of 
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1000 km s -1 . In other words, a high increase in 
the velocity dispersion of the stars toward the dy- 
namical center is revealed. Furthermore, a large 
and coherent counter-rotation, especially of the 
early-type stars, is shown, supporting their origin 
in a well-defined epoch of star formation. Obser- 
vations of stellar winds nearby Sgr A* give a mass 
accretion rate of dM/dt = 6 x 1O~ 6 M yr" 1 (Gen- 
zel et al. 1996). Hence, the dark mass must have 
a density ~ 10 9 Mq pc -3 or greater, and a mass- 
to-luminosity ratio of at least 100M Q /i Q . The 
bottom line is that the central dark mass seems to 
be a single object, and that it is statistically very 
significant (<~ 6 — 8cr). 

2.2. Domain of the black hole? 

Such a large density contrast excludes that the 
dark mass could be a cluster of almost 2 x 10 6 
neutron stars or white dwarfs. Detailed calcu- 
lations of evaporation and collision mechanisms 
give maximal lifetimes of the order of 10 s years, 
much shorter than the estimated age of the Galaxy 
(Sanders 1992; Maoz 1995). 

As a first conclusion, several authors state that 
in the Galactic Center there is either a single su- 
permassive black hole or a very compact cluster 
of stellar size black holes (Genzel et al. 1996). We 
shall come back to this paradigm through the rest 
of the paper (particularly in Section III). 

Dynamical evidence for central dark objects has 
been published for 17 galaxies, like M87 (Ford 
et al. 1994; Macchetto et al. 1997), or NGC4258 
(Greenhill et al. 1995), but proofs that they re- 
ally are black holes requires measurement of rel- 
ativistic velocities near the Schwarzschild radius, 
r s ~ 2M./(10 8 M Q ) AU (Kormcndy and Rich- 
stone 1995), or seeing the properties of the ac- 
cretion disk. Due to the above mentioned mass 
accretion rate, if Sgr A* is a supermassive black 
hole, its luminosity should be more than 10 40 erg 
s _1 , provided the radiative efficiency is about 10%. 
On the contrary, observations give a bolometric 
luminosity less than 10 37 erg s _1 , already taking 
into account the luminosity extinction due to in- 
terstellar gas and dust. This discrepancy is the 
so-called "blackness problem" which has led to the 
notion of a "black hole on starvation". Standard 
dynamics and thermodynamics of the spherical ac- 
cretion onto a black hole must be modified in or- 
der to obtain successful models (Falcke and Hein- 



rich 1994). Recent observations, we recall, probe 
the gravitational potential at a radius larger than 
4 x 10 4 Schwarzschild radii of a black hole of mass 
2.6 x 10 6 M Q (Ghez 1998). 

2.3. Domain of the massive neutrino? 

An alternative model for the supermassive com- 
pact object in the center of our Galaxy has been 
recently proposed by Tsiklauri and Viollier (1998). 
The main ingredient of the proposal is that the 
dark matter at the center of the galaxy is non- 
baryonic (but in any case fermions, e.g. massive 
neutrinos or gravitinos), interacting gravitation- 
ally to form supermassive balls in which the de- 
generacy pressure balances their self-gravity. Such 
neutrino balls could have formed in early epochs, 
during a first-order gravitational phase transition 
and their dynamics could be reconciled, with some 
adjustments, to the Standard Model of Cosmology. 

Several experiments are today running to 
search for neutrino oscillations. LSND (White 
1998) found evidence for oscillations in the v e — 
channel for pion decay at rest and in flight. On 
the contrary KARMEN (Zeitnitz 1998) seems to 
be in contradiction with LSND evidence. CHO- 
RUS and NOMAD at CERN are just finishing the 
phase of 1994-1995 data analysis. It is very likely 
that exact predictions for and — v T oscillations 
will be available in next years. Thanks to the fact 
that it is possible to give correct values for the 
masses to the quarks up, charm, and top, it is 
possible to infer reasonable values of mass for v e , 
Vp, and v T . In order to explain the characteristics 
of Sgr A*, we need fermions whose masses range 
between 10 and 25 keV which, cosmologically, fall 
into the category of warm dark matter. It is in- 
teresting to note that a good estimated value for 
the mass of r-neutrino is 

m„ T = ro„„ ( — ) ~ 14.4 keV. (1) 

Choosing fermions like neutrinos or gravitinos in 
this mass range allows for the formation of su- 
permassive degenerate objects from 10 6 M Q to 
1O 9 M . 

The theory of heavy neutrino condensates, 
bound by gravity, can be easily sketched consider- 
ing a Thomas-Fermi model for fermions (Viollier 
et al. 1992). We can set the Fermi energy equal 
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to the gravitational potential which binds the sys- 
tem, and see that the number density is a function 
of the gravitational potential. Such a gravitational 
potential will obey a Poisson equation where neu- 
trinos (and anti-neutrinos) are the source term. 
This equation is valid everywhere except at the 
origin. By a little algebra, the Poisson equation 
reduces itself to the radial Lane-Emden differ- 
ential equation with polytropic index n = 3/2, 
which is equivalent to the Thomas-Fermi differen- 
tial equation of atomic physics, except for a minus 
sign that is due to the gravitational attraction of 
the neutrinos. The Newtonian potential, which is 
the solution of such an equation, is <&(r) ~ r~ 4 . 
Considering a standard accretion disk, if Sgr A* 
is a neutrino star with radius R = 30.3 Ids (<~ 10 5 
Schwarzschild radii), mass Mac = 2.6 x 1O 6 M , 
and luminosity Lqc ~ 10 37 erg sec -1 , it should 
consist of neutrinos with masses m > 12.0 keV 
for g v — 4, or m > 14.3 keV for g v = 2 (Bilic 
1998; Viollier et al. 1992). It is specially ap- 
pealing that with the same mass for neutrinos, 
several galactic centers can be modeled. For in- 
stance, similar results hold also for the dark object 
(M~3x 10 9 M Q ) inside the center of M87. Due 
to the Thomas-Fermi theory, the model fails at 
the origin: we have to consider the effect of the 
surrounding baryonic matter which, in some sense, 
has to stabilize the neutrino condensate. In fact, 
the solution $(r) ~ r~ 4 is clearly unbounded from 
below. 

The model proceeds assuming a thermodynam- 
ical phase where a constant neutrino number den- 
sity can be taken into consideration. This is quite 
natural for a Fermi gas at temperature T = 0. 
The Poisson equation can be recast in a Lane- 
Emden form with polytropic index n = 0. The 
solution of such an equation is $(r) ~ r 2 , which 
is clearly bounded (Capozziello and Iovane 1999). 
For T ^ 0, we get a solution of the form $(r) ~ 
r~ 4 . Matching these two results it is possible to 
confine the neutrino ball. On the other hand, a 
similar result is recovered using the Newton The- 
orem for a spherically symmetric distribution of 
matter of radius R (Binney and Tremaine 1987). 
In that case, the potential goes quadratic inside 
the sphere while it goes as $(r) ~ r -1 match- 
ing on the boundary. In our case, the situation 
is similar assuming the matching with a steeper 
potential. 



If such a neutrino condensate exists in the cen- 
ter of Galaxy, it could act as a spherical thick lens 
(a magnifying glass) for the stars behind it, so that 
their apparent velocities will be larger than in real- 
ity. In other words, depending on the line of sight, 
it should be possible to correct the projected ve- 
locities by a gravitational lensing contribution, so 
trying to explain the bimodal distribution (early 
and late type stars) actually observed (Ghez 1998; 
Genzel et al. 1996). Since the astrophysical fea- 
tures of the object in Sgr A* are quite well known, 
accurate observations by lensing could contribute 
to the exact determination of particle constituents. 
A detailed model and comparison with the data 
was presented by some of us (Capozziello and Io- 
vane 1999). 

3. Observational status 

3.1. Brief review on dynamical data and 
what are they implying 

An important information on the central ob- 
jects of galaxies (particularly in active galactic nu- 
clei) is the short timescale of variability. This has 
the significance of putting an upper bound -known 
as causality constraint- on the size of the emitting 
region: If a system of size L suddenly increases its 
emissivity at all points, the temporal width with 
which we receive it is L/c and thus, a source can 
not fluctuate in a way that involves its entire vol- 
ume in timescales shorter than this (unless c is not 
a limiting velocity). This finally yields a corre- 
sponding maximum lenghscale, typically between 
10~ 4 to 10 pc (Krolik 1999). Autocorrelation in 
the emission argue against the existence of a clus- 
ter of objects (unless only one member dominates 
the emissivity), and the small region in which the 
cluster should be allows two body encounters to 
be very common and produces the lost of the sta- 
bility of complex clusters. These facts can be used 
to conclude that a single massive object must be 
in the center of most galaxies. What observations 
show in this case is that the size of the emitting 
regions are very small. This does not directly im- 
plicate black holes as such. 

The way in which we expect to detect the influ- 
ence of a very massive object is through its gravity. 
The "sphere of influence" of a large mass is defined 
by the distance at which its potential significantly 
affects the orbital motions of stars and gas, and is 
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jiven by 



R* = GM/a 2 ~ 4M 7 aJ 00 pc, 



(2) 



where the central mass is normalized to 10 7 Mq 
and ct^^oo is a rms orbital speed in units of 100 
km s _1 . Thus, even very large masses have a small 
sphere of influence. Within this sphere, the ex- 
pected response to a large nuclear mass can be 
divided in two groups. Firstly, the response of in- 
terstellar gas can hardly be anything different from 
an isotropic motion in the center of mass frame. 
Random speeds are often thought to be much 
smaller than this overall speed, which far from the 
sphere of influence is just given by (GM /r) 1 / 2 . If 
the energy produced by random motions is radi- 
ated away, the gas behaves as a whole and tend to 
flatten itself, conserving the angular momentum. 
Observations of the rotational velocity of gas as a 
function of the radius can thus provide a measure 
of the total mass at the center. Then, what ob- 
servation searches is the existence of a Keplerian 
potential signature in the flattened gas velocity dis- 
tribution. Examples of this are the radio galaxies 
M87 (Ford et al. 1994), NGC 4258 (Miyoshi et al. 
1995), and many others. Then, any massive object 
producing a velocity distribution with a Keplerian 
decrease would be allowed by observation. 

Secondly, we consider the response of the stars. 
In this case, peculiar velocities are larger or com- 
parable to the bulk streaming and it is harder to 
actually differentiate the stellar response to a nu- 
clear mass from the observable properties of a pure 
stellar potential. However, stellar motions, con- 
trary to that of gas, is not affected by other forces 
(as those produced by magnetic fields), and are 
better tracers of mass distributions. What obser- 
vations show in this case are the stellar density 
and the velocity distribution. Use of the collision- 
less Boltzmann equation allows one to get Jean's 
relation (Binney and Tremaine 1987; Krolik 1999): 



GM(r) 



d\nn{r) d In a 2 
d In r dh\r 



+ 2- 



(3) 

here, n(r) is the spherically symmetric distribu- 
tion of stars, V rot is the mean rotational speed, 
Pr,e,<j> are velocity dispersions and of = a 2 + o|. 
Measurements of V and random velocities deter- 
mine that M(r) decline inwards until a critical ra- 
dius, where M becomes constant. What really 



happens, however, if that even if this region is not 
scrutinized, the mass to light ratio becomes suf- 
ficiently large to suggest the presence of a large 
dark mass. A combination of this with other mea- 
surements strongly suggest that this dark mass is 
a single object. 

To summarize, usual measurements point to 
inform us that there exist a single supermassive 
object in the center of some galaxies but do not 
definitively state its nature. As Kormcndy and 
Richstone (1995) have stated, the black hole sce- 
nario has become our paradigm. But suggestions 
that the dark objects are black holes are based 
only on indirect astrophysical arguments, and sur- 
prises are possible on the way to the center. 

3.2. The Galaxy 

We follow Genzel et al. (1996) and parameter- 
ize the stellar density distribution as, 



n(r) 



S 



Ro 1 + (R/R ) a ' 



(4) 



Note that Ro is related to the core radius through 
-Rcore = b(a)Ro. Genzel et al. found that the 
best fit parameters for the observed stellar cluster 
are a central density of 4 x 10 6 -M©pc~ 3 , a core 
radius of 0.38 pc, and a value a = 1.8 (b(a) = 
2.19). With this distribution, they found that a 
dark mass of about 2.5 xl0 6 M Q was needed to fit 
the observational data. 

The cluster distribution, and the cluster plus 
a black hole constant mass is shown in Fig. 1. 
Black boxes represent Genzel et al. 1996 (Table 
10) and Eckart and Genzel 1997 (Fig. 5) data. A 
dashed line stands for the stellar cluster contribu- 
tion, while a dot-dashed line represent an enclosed 
point-like black hole mass. The solid line in the 
right half of the figure stands for the mass dis- 
tribution both for a black hole and a boson star 
plus the stellar cluster. The mass dependence we 
are plotting for the boson star is obtained in Sec- 
tion IV, and represents the mass distribution of a 
mini-boson star (A = 0, m[GeV] = 2.81 x 10~ 26 ) 
with dimensionless central density o(0) equal to 
0.1. Other boson star configurations, with appro- 
priate choice of the boson field mass to, yield to the 
same results. Note the break of more than three 
orders of magnitude in the x-axis, this is caused 
because the boson star distribution, further out of 
the equivalent Schwarzschild radius, behaves as a 
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black hole. It begins to differ from the black hole 
case at radius more than three orders of magnitude 
less than the innermost data point, that is why the 
x-axis has to have a break. From the mass distri- 
bution, a boson star in the center of the galaxy is 
virtually indistinguishable from a black hole. 

Tsiklauri and Viollicr (1998) have shown that 
the same observational data can also be fitted us- 
ing an extended neutrino ball. In that case, dif- 
ferences begin to be noticed just around the in- 
nermost data point. It is then hard to determine 
whether the central object is a black hole, a neu- 
trino ball, or a boson star based only on dynamical 
data now at hand, and other problems concerning 
the accretion disk have to considered (see below). 
Even harder is the situation for deciding -using 
only this kind of data- if the supermassive object 
is a boson star instead of a black hole: as a boson 
star is a relativistic object, the decay of the en- 
closed mass curve happens close the center. This, 
however, will provide an equivalent picture than a 
black hole for disruption and accretion processes; 
we shall comment on it in the next sections. 

To apply Eq. (3) to the observational data 
we have to convert intrinsic velocity dispersions 
(oy(r)) and volume densities (n(r)) to projected 
ones (Binney and Tremaine 1987), these are the 
ones we observe. We shall also consider that 
(5 = 1 — <jg/(Tr, the anisotropy parameter, is equal 
to 0, and we are assuming implicitly that <jg — a^. 
We take into account the following Abel integrals, 

S(p) denotes surface density, <r r (p) is the pro- 
jected velocity dispersion, and p is the projected 
distance. We adopt Genzel et al.'s (1996) and 
Tsiklauri and Viollicr's (1998) parameterization 
for oy(r), 

<7 r (r)=a(oo) 2 +<7(2") 2 ^J . (7) 

What one usually does is to numerically inte- 
grate Eqs. (5,6) and fit the observational data 
(a r (p) vs. p). To do so one also has to assume a 
density dependence for the cluster (as in Eq. (4)), 



and the dark mass. With a point-like dark mass 
of 2.5 xl0 6 M Q , the parameters of the fit result to 
be cr(oo) = 55 km s" 1 , a (2") = 350 km s"\ and 
j3 = 0.95. 

If one now changes the central black hole for 
a neutrino ball, or a boson star, one has to con- 
sider the particular density dependence for these 
objects. In this way, n(r) = n stcllar clustor (r) + 

T^dark mass 

(r). We obtained riboson star (r) as M(r)/(4/37rr 3 ), 
where M(r) is the fitted mass dependence of the 
boson star as explained in the next section. It 
is now useful to consider that the fitting of a r (p) 
is made taking into account observational data 
points in regions where the boson star generated 
space-time is practically indistinguishable from a 
black hole. Then, we may expect that the actual 
parameters, er(oo), er(2"), and /?, will be very close 
to those obtained for the black hole. For our pur- 
poses, it is enough to take the same a r (r) as in 
the black hole case, and compute ay boson star (p) 
using Eqs. (5,6), with the adequate total n(r). 

In Fig. 2 we show the observational data 
of Eckart and Genzel (1997) (filled black boxes) 
and Genzel et al. (1996) (hollow circles) super- 
imposed with the curve for <J r (p) that we ob- 
tain with a mini-boson star (A = 0, <r(0) = 0.1, 
m[GcV] = 2.81 x 10~ 26 ) in the galactic center. 
Other boson star configurations, with appropri- 
ate choice of the boson field mass m, yield to 
the same results. For this configuration, we used, 
as will be explained in Section IV, a mass distri- 
bution given by a Boltzmann-like equation with 
A 2 = 2.5 x 1O 6 M , Ax = -0.237 x 1O 6 M , 
Ra = 1.19210" 6 pc, and AR = 4.16310- 6 pc. In 
the range plotted, and in which data is available, 
the differences between boson and black hole theo- 
retical curves is undetectable. They only begin to 
deviate from each other at p ~ 10 pc, well be- 
yond the last observational data point. Even the 
deviation in such a region is as slight as 1 km s _1 , 
and it only becomes more pronounced when p val- 
ues are closer to the center. However, we should 
recall that it has no sense to go to such extreme 
values of p: the stars will be disrupted by tidal 
forces (see the discussion in Section VI) in those 
regions. 
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4. Boson stars 

4.1. Basic concepts and configurations 

Let us study the Lagrangian density of a mas- 
sive complex self-gravitating scalar field (taking 

h= c= 1), 



C = 



in 



(8) 

where R is the scalar of curvature, \g\ the modulus 
of the determinant of the metric g M „, and ip is a 
complex scalar field with potential U. Using this 
Lagrangian as the matter sector of the theory, we 
get the standard field equations, 



R u 



1 

dU 



rn 



Pi 



0, 



(9) 
(10) 



where the stress energy tensor is given by, 
1 



g a P{d a r){d^)-U{\^\]l) 



and 



□ = d» 



9 I <Td v I 



(12) 



is the covariant d'Alembcrtian. Because of the 
fact that the potential is a function of the square 
of the modulus of the field, we obtain a global J7(l) 
symmetry. This symmetry, as we shall later dis- 
cuss, is related with the conserved number of par- 
ticles. The particular form of the potential is what 
makes the difference between mini-boson, boson, 
and soliton stars. Conventionally, when the po- 
tential is given by 



U = m 2 \yj\ 2 + ^ 



(13) 



where m is the scalar mass and A a dimensionless 
constant measuring the self-interaction strength, 
mini-boson stars are those spherically symmetric 
equilibrium configurations with A = 0. Boson 
stars, on the contrary, have a non-null value of A. 
The previous potential with A =/= was introduced 
by Colpi et al. (1986), who numerically found that 
the masses and radius of the configurations were 



deeply enlarged in comparison to the mini-boson 
case. 

Soliton (also called non-topological soliton) 
stars arc different in the sense that, apart from 
the requirement that the Lagrangian must be in- 
variant under a global U(l) transformation, it 
is required that -in the absence of gravity- the 
theory must have non-topological solutions; i.e. 
solutions with a finite mass, confined to a finite 
region of space, and non-dispersive. An example 
of this kind of potentials is the one introduced by 
Lee and his coworkers (Friedberg et al. 1987), 



£/ = m 2 H 2 



1 - 



IVf 



(14) 



where $ is a constant. In general, boson stars 
accomplish the requirement of invariance under a 
U(l) global transformation but not the solitonic 
second requirement. To fulfill it, it is necessary 
that the potential contains attractive terms. This 
is why the coefficient of (4>*%li) 2 of the Lee form has 
a negative sign. Finally, when |^| — > oo, U must 
be positive, which leads, minimally, to a sixth or- 
der function of ip for the self-interaction. It is usu- 
ally assumed, because of the range of masses and 
radius for soliton stars in equilibrium, that they 
are huge and heavy objects, although this finally 
depends on the choice of the different parameters. 

We shall now briefly explain how these con- 
figurations can be obtained (Kaup 1968; Ruffini 
and Bonazzola 1969; Lee and Pang 1992; Liddle 
and Madsen 1992; Mielke and Schunck 1998). We 
adopt a spherically symmetric line clement 

ds 2 = e u{r) dt 2 - e^dr 2 - r 2 {d$ 2 + sin 2 i? dip 2 ) , 

(15) 

with a scalar field time-dependence ansatz consis- 
tent with this metric: 



il>(r, t) = a(r)e 



-iujt 



(16) 



where u is the (eigen-)frcquency. This form of the 
field ensures us to be working in the configurations 
of minimal energy (Friedberg et al. 1987). 

The non-vanishing components of the energy- 
momentum tensor are 

T ° = p = l -[uJ 2 a 2 {r)e-» + a' 2 (r)e^ + U] , (17) 
Ti 1 =Vr = \[u 2 o 2 {r)e- v + o l2 {r)e-» - U] , (18) 
T 2 2 = T 3 3 =p ± = - X -[u 2 o 2 {r)e- v - a' 2 (r)e^ - U] (19) 
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where ' = d/dr. One interesting characteristic 
of this system is that the pressure is anisotropic; 
thus, there are two equations of state p r = p — U 
and p± = p — U — cr' 2 (r)e _Ai . The non- vanishing 
independent components of the Einstein equation 
are 



u' + p' 



8tt 



m 



(p + p r )re^ 



(20) 



pi 



^pre* - -(e" - 1) . (21) 
m.pj r 



Finally, the scalar field equation is 



„ , , v> - pi 2 



+ - I a' + e^~ v w 2 a - e^^—a = . 

da z 

(22) 

To do numerical computations and order of mag- 
nitude estimates, it is useful to have a new set of 
dimensionlcss variables. We adopt here 



x = mr, 



(23) 



for the radial distance, we redefine the radial part 
of the boson field as 



a = V47T cr/mpi, 



and introduce 



A = \mpj4Trm 2 , fi = — . 

m 



(24) 



(25) 



In order to obtain solutions which are regular at 
the origin, we must impose the following boundary 
conditions cr'(0) = and p(0) — 0. These solu- 
tions have two fundamental parameters: the self- 
interaction and the central density (represented by 
the value of the scalar field at the center of the 
star). The mass of the scalar field fixes the scale 
of the problem. Boundary conditions representing 
asymptotic flatness must be applied upon the met- 
ric potentials, these determine -what is actually 
accomplished via a numerical shooting method- 
the initial value of v = f(0). Then, having de- 
fined the value of the self interaction, or alterna- 
tively, the form of the soliton potential, the equi- 
librium configurations are parameterized by the 
central value of the boson field. As this central 
value increases, so does the mass and radius of the 
the star. This happens until a maximum value 
is reached in which the star looses its stability 
and disperses away (the binding energy being pos- 
itive). Up to this value of cr , catastrophe theory 



can be used to show that these equilibrium config- 
urations are stable (Kusmartsev et al. 1991). As 
an example of boson star configurations, we show 
in Fig. 3, the mass and number of particles (see 
below) for a A = 10 Colpi et al.'s potential, and 
in Fig. 4, the stability analysis. 

When A ^> 1 , as in some of the cases we explore 
in the next sections, we must follow an alternative 
adimensionalization (Colpi et al. 1986). For large 
A, we shift to the following set of variables, 



<r. = ok 1 ' 2 , 

= xk- 1 ' 2 , 
M* = MA" 1/2 . 
Here, is defined by 



(26) 

(27) 
(28) 

(29) 



which corresponds to the Schwarzschild mass (see 
below). Ignoring terms C(A _1 ), the scalar wave 
equation is solved algebraically to yield, 



(30) 



and up to the same accuracy, the field equations 
are 

^l = l x l(3n 2 e-» + l)(n 2 e-»-l), (31) 
ax* 4 

-P— ~ {l-e-») = \{tfe-»-lf. (32) 

The system now depends only on one free param- 
eter Q 2 e~"^ . Numerical solutions show that the 
maximum mass corresponding to a stable star is 
given by M max - 0.22A 1 / 2 m 2 31 /m. 

4.2. Masses estimates 

The invariance of the Lagrangian density un- 
der a global phase transformation ip — > ipe^' 1 " 9 of 
the complex scalar field gives (via the Noether's 
theorem) a locally conserved current d^j^ = 0, 
and a conserved charge (number of particles). We 
need to study the number of particles because it 
is essential to determine whether the configura- 
tions are stable or not. A necessary requirement 
towards the stability of the configurations is a neg- 
ative binding energy (BE = M — mN), i.e. the 
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star must be energetically more favorable than a 
group of unbound particles of equal mass. From 
the Noether theorem, the current is given by 

f = \VUT\g^[rd^ - i>d v p] • (33) 

and the number of particles is 

N := J fd 3 x. (34) 

For the total gravitational mass of localized solu- 
tions, we may use Tolman's expression (Tolman 
1934), or equivalently, the Schwarzschild mass: 



M 



= J(2T °-Tf)^\J}d 3 



(35) 



In Fig. 5, we show the mass of a boson star 
with A = as a function of x. We have fitted this 
curve with a Boltzmann-like function 



M m (x) = A 2 



1 + fj 1 ^ 10 )/ 41 ' 



(36) 



which is reliable except in regions very near the 
center, where M& t {x) becomes slightly negative. 
The x 2 -paramcter of the fitting is around 10~ 5 , 
and values for A\,A<i, and Ax are given in the fig- 
ure. It is this formula what we have used in Fig. 
2 to analytically get the a r (p) dependence of the 
boson star (in the range we use the approximation, 
the actual mass and the fitting differ negligibly). 
Note also what the fitting is physically telling us: 
it represents a black hole of mass A\ plus an in- 
ner exponentially decreasing correction. This is 
the first time that such a Boltzmann-like fitting is 
done, and we think it could be usefully applied in 
other analytical computations. Models with dif- 
ferent A and a can be equally well fitted. 

Since boson stars are prevented from gravita- 
tional collapse by the Heisenberg uncertainty prin- 
ciple, we may make some straightforward mass 
estimates (Mielke and Schunck 1998): For a bo- 
son to be confined within the star of radius Ro, 
the Compton wavelength has to satisfy A^, = 
(2irh/mc) < 2R . In addition, the star radius 
must be of the order of the last stable Kepler or- 
bit 3i?s around a black hole of Schwarzschild ra- 
dius i?s := 2GM. In the case of a mini-boson 
star of effective radius Ro = (ir/2) 2 Rs close to its 
Schwarzschild radius one obtains the estimate 

M crit = (2/ 7 r)m| 1 /m > 0.633 M^/m . (37) 



The exact value in the second expression was 
found only numerically. 

For a mass of m = 30 GcV, one can estimate the 
total mass of this mini-boson star to be M ~ 10 10 
kg and its radius i? — 10~ 17 m, amounting a 
density 10 48 times that of a neutron star. In the 
case of a boson star (A ^ 0), since |V>| ~ mpi/\/87r 
inside the boson star (Colpi et al. 1986), the energy 
density is 



p ~ rr?m% x (1 + A/8) 



(38) 



Equivalently, we may think that this corresponds 
to a star formed from non-interacting bosons with 
re-scaled mass m —> m/ ^/T+ A/8 and conse- 
quently, the maximal mass scales with the cou- 
pling constant A as, 



- in 



M crit ~ -Vl + A/8 

ir m 



pi 



(39) 



There is a large range of values of mass and radius 
that can be covered by a boson star, for different 
values of A and ct - F° r instance, if m is of the 
order of the proton mass and A ~ 1, this is in the 
range of the Chandrasekhar limiting mass Mch '■= 
M^/m 2 ~ 1.5M . 

Larger than these estimates is the range of 
masses that a non-topological soliton star pro- 
duces, this is because the power law depen- 
dence on the Planck mass is even higher: <~ 
10 _2 (TOpj/m$Q). These configurations are static 
and stable with respect to radial perturbations. 
It is by no means true, however, that these are 
the only stable equilibrium configurations that 
one can form with scalar fields. Many extensions 
of this formalism can be found. It is even pos- 
sible to see that boson stars are a useful setting 
where to study gravitation theory in itself (Torres 
1997; Torres et al. 1998c,d; Whinnett and Tor- 
res 1999). Most importantly for our hypothesis 
is that rotating stable relativistic boson stars can 
also be found with masses and radii comparable in 
magnitude to their static counterparts (Schunck 
and Mielke 1996; Yoshida and Eriguchi 1997). In 
astrophysical settings it is usual to expect some 
induced rotation of stellar objects, and it is im- 
portant that these rotation may not destabilize 
the structure. Other interesting generalization is 
that of electrically charged boson stars introduced 
by Jetzer and van der Bij (1989). Although it 
is usually assumed that selective accretion will 
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quickly discharge any astro-physical object, some 
recent results by Punsly (1998) suggest this may 
not always be the case. 

From these simple considerations, we can set 
a difference between boson and fermion conden- 
sations (e.g. neutrino condensation): in the sec- 
ond case we can have an extended object (in the 
case of Sgr A*, its size can be of about <~ 301d) 
and it can be quite diluted. In the case of a bo- 
son star, the object is "strongly" relativistic, ex- 
tremely compact, and its size is comparable with 
its Schwarzschild radius. We shall discuss the con- 
sequences of this point widely in the next section. 

4.3. Effective potentials 

The motion of test particles can be obtained 
from the Euler-Lagrange equations taking into ac- 
count the conserved canonical momenta (Shapiro 
and Teukolsky 1983). In a general spherically sym- 
metric potential, the invariant magnitude squared 
of the four velocity (it 2 = g^i^x"), which is 1 for 
particles with non-zero rest mass and for mass- 
less particles, yields to 



r 2 = 



9tt 



9<t»t> 



(40) 



where I and are constants of motion given by 
— 9<P<P sin 2 9 <fi and gut (angular momentum and en- 
ergy at infinity, both per unit mass) respectively, 
and a dot stands for derivation with respect to 
an affine parameter. This equation can be trans- 
formed to 

\f- 2 + V cS = \El e-^\ (41) 

where V c s is an effective potential. This name 
comes from the fact that in the Schwarzschild so- 



lution 



1 and thus, the previous equation 



can be understood as a classical trajectory of a 
particle of energy E 2 /2 moving in a central po- 
tential. This is not so in a more general spherically 
symmetric case, like these non-baryonic stars. In 
the boson star case, for instance, typical met- 
ric potentials are shown in Fig. 6; note that 
for them e~^~ v ^ 1 We can see, however, that 
e -v.-v < e -»(0)-v(0) = Cj where c is a cons t a nt, 

and then, usual classical trajectories can be looked 
at, in the sense that we may construct an equation 
of the form r 2 /2 + V cS < l/2E% c e-^-» ( - \ and 
it will be always satisfied. 



The effective potential that massive or massless 
particles would feel must be very different from the 
black hole case. In the case of a massless particle, 
the effective potential is given by 



while for a massive test particle it is 
V cS 



1 „ / I 2 
■ze-» 1 + 



(42) 



(43) 



We show both, boson (for the case of A = 0) and 
black hole potentials in Figs. 7 and 8. The mass 
of the central object is fixed to be the same in 
both cases (see captions) and the curves represent 
a fourth order Runge-Kutta numerical integration 
of the equations we previously derived with A = 
and c(0) = 0.1. In the case of massive test parti- 
cles we use l 2 /M 2 = 0, 12 and 15; thus explicitly 
showing the change in the behavior of V e s for the 
black hole case. As / increases, this shape changes 
from a monotonous rising curve to one that has 
a maximum and a minimum before reaching its 
asymptotic limit. These extrema disappear for 
l 2 /M 2 < 12. In the case of boson stars, however, 
if I 7^ we have a divergence in the center at r = 
and only one extremum -a minimum-, which oc- 
curs at rising values ofr/M as I grows. The curve 
V e ff for I = -radially moving objects- is not diver- 
gent, and we can see that the particles may reach 
the center of the star with a non-null velocity, 
given by l/2r 2 = l^E^e'^-^ - V eS (0) > 0, 
and will then traverse the star unaffected. 

For massless particles, differences are also noto- 
rious: a black hole produces a negative divergence 
and a boson star a positive one. Radial motion 
of massless particles is insensitive to V e s, being 
this equal to zero, as in the Schwarzschild case. 
In both cases, for massive and massless particles, 
outside the boson star the potential mimics the 
Schwarzschild one. 

4.4. Particle orbits 

In the case of a black hole, we can use the New- 
tonian analogy. Orbits can be of three types: if 
the energy is bigger than the effective potential at 
all points, particles are captured. If the energy is 
such that the energy equals the effective potential 
just once, particles describe an unbound orbit, and 
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the point of equality is known as turning point. If 
there are two such points, orbits are bound around 
the black hole. Orbits in which the energy equals 
the potential in a minimum of the latter are cir- 
cular and stable (r = 0, < 0). 

We can not, because of the different relation- 
ship between the metric potentials we commented 
above, do the same analysis using V e g for the bo- 
son star. We can, however, note that in most cases 
the total equation for the derivative of r will be 
modified in a trivial way: If we look at the cases 
where the effective potential has a divergence at 
the center (and because the metric coefficients do 
not diverge), there is no other possibility for the 
particles more than found a turning point. Or- 
bits can then be bound or unbound depending on 
the energy, but we can always find a place where 
r = and then it has to reverse its sign. In the 
particular cases in which the equality happens at 
the minimum of the potential we have again sta- 
ble circular orbits. Only in the case where I = 
particles can traverse the scalar star unaffected. 
For m 7^ there is still the possibility of finding a 
turning point, if the energy is low enough. How- 
ever, if the particle is freely falling from infinity, 
with Eoo — 1, all energy is purely rest mass, it will 
radially traverse the star, as would do a photon. 

We conclude that all orbits are not of the cap- 
ture type. They can be circular, or unbound, and 
they all have at least one turning point. This 
helps to explain why a non-baryonic object will 
not develop a singularity while still being a rel- 
ativistic object (comparable effective potentials, 
equivalently, comparable relativisticity coefficient 
GM(r)/r.). 

5. Galactic parameters and mass 

One interesting fact, which seems to be not 
referred before, is the point that for all these 
scalar stars, their radius is always related with 
the mass in the same way: R > Mrripf. This 
is indeed the statement -contrary to what it is 
usually assumed- that not all interesting astro- 
physical ranges of mass and radius can be mod- 
eled with scalar fields. In the scalar star models, 
from the given central mass, the radius we obtain 
for the star is comparable to that of the horizon, 
R ~ to p1 2 x 2.61 x 10 6 M Q - 3.9 x 10 n cm. 

The question now is for which values of the pa- 



rameters we can obtain a scalar object of such a 
huge mass. For the case of mini-boson star we 
need an extremely light boson: 



m[GcV] = 1.33 x 10 



-25 



Af(oo) 



BH 



(44) 



Given a central density cr(0), M(oo) stands for the 
dimcnsionlcss value of the boson star mass as seen 
by an observer at infinity. Mbh is the value of 
the black hole mass (in millions of solar masses), 
obtained by fitting observational data. Then, wc 
are requiring that the total mass of the boson star 
equals that of the black hole. For instance, in Fig. 
1, we have taken Eckart and Genzel's (1997) and 
Genzel et al.'s (1996) data, and fit them with a 
mini-boson star with ct(0) = 0.1, which yields to a 
boson mass given by m[GeV] = 2.81 x 10~ 26 ; the 
total mass of the star (without the cluster contri- 
bution) is 2.5 xlO 6 M Q . 

In the case of boson stars, and using the critical 
mass dependence, oc ^/Xm^/m 2 , the requirement 
of a 2.5 million mass star yields to the following 
constraint, 



m[GeV] = 7.9 x 10~ 4 ( A 



1/4 



(45) 



It is possible to fulfill the previous relationship, 
for instance, with a more heavy boson of about 
1 MeV and A ~ 1. A plot of this relation -for 
some values of A- is shown in Fig. 9. Note that in 
this case, the value of the dimcnsionlcss parameter 
A is huge, and special numerical procedures, as 
explained above, must be used to obtain solutions. 
The characteristics of these solutions have proven 
to be totally similar to those with A = which 
were used in Fig. 1, just the adimensionalization 
differs. 

Finally, in the case of a non-topological soliton 
we obtain the following constraint, 



m[GeV] 



7.6 x 10 12 



(46) 



For the usually assumed case, in which the or- 
der parameter <&o is of equal value than the bo- 
son mass, we need very heavy bosons of unit mass 
m = 1.2 x 10 4 GeV. Other possible pairs are shown 
in Fig. 10. 



11 



5.1. Boson candidates? 

Based only on the constraints imposed by the 
mass-radius relationship valid for the scalar stars 
analyzed, we may conclude that: 

1. if the boson mass is comparable to the 
expected Higgs mass (hundreds of GeV), 
then the Center of Galaxy could be a non- 
topological soliton star; 

2. an intermediate mass boson could produce 
a super-heavy object in the form of a boson 
star; 

3. for a mini-boson star to be used as central 
objects for galaxies it is needed the existence 
of an ultra-light boson. 

These conclusions should be considered as order 
of magnitude estimations. Several reasons force 
us to make this warning. Firstly, we are just con- 
sidering static and uncharged stars, this is just a 
model (the simplest), but more complicated ones 
can modify the actual constraints. Secondly, we 
do not know the exact form of the self-interaction, 
or in the case of non-topological stars, the value 
of $o. For instance, consider the Higgs mass. 
In the electro- weak theory a Higgs boson doublet 
(<j> + ,$°) and its anti-doublet ($~,$°) are neces- 
sary ingredients in order to generate masses for the 
W ± and Z° gauge vector bosons. Calculations of 
two-loop electro-weak effects have lead to an in- 
direct determination of the Higgs mass (Gambino 
1998). For a top quark mass of M t = 173.8 ± 5 
GeV, the Higgs mass is mj, = 104tg GeV. How- 
ever, experimental constraints are weak. Fermi- 
lab's tevatron (Han and Zhang 1999) has a mass 
range of 135 < < 186 GeV, and together with 
LHC at Cern, they could decide if these Higgs par- 
ticles (in the given range) exist in nature. In- 
teresting is to note that as a free particle, the 
Higgs boson is unstable with respect to the de- 
cays h —> W + + W~ and h -> Z° + Z°. However, 
as remarked by Mielke and Schunck (2000), in a 
compact object, and in full analogy with neutron 
stars -where there is equilibrium of (3 and inverse j3 
decay- these decay channels are expected to be in 
equilibrium with the inverse process Z° + Z a — ► h. 

We should also mention the possible dilatons 
appearing in low energy unified theories, where 
the tensor field g^ v of gravity is accompanied 



by one or several scalar fields. In string effec- 
tive super-gravity (Ferrar et al. 1994), for in- 
stance, the mass of the dilaton can be related 
to the super-symmetry breaking scale m susy by 
m v ~ 10~ 3 (m susy / TeV) 2 eV. Finally, a scalar 
with a long history as a dark matter candidate 
is the axion, which has an expected light mass 
m a = 7.4 x (10 7 GeV// CT ) eV > 10" 11 eV with de- 
cay constant j a close to the inverse Planck time. 
Goldstone bosons have also inferred mass in the 
range of eV and less, m g < 0.06 — 0.3 eV (Umeda 
et al. 1998). 

If boson stars really exist, they could be the 
remnants of first-order gravitational phase tran- 
sitions and their mass should be ruled by the 
epoch when bosons decoupled from the cosmolog- 
ical background. The Higgs particle, besides its 
leading role in inflationary theories, should be the 
best and natural candidate as constituent of a bo- 
son condensation if the phase transition occurred 
in early epochs. A boson condensation should be 
considered as a sort of topological defect relic. In 
this case, as we have seen, Sgr A* could be a soli- 
ton star. If soft phase-transitions took place dur- 
ing cosmological evolution (e.g. soft inflationary 
events), the leading particles could have been in- 
termediate mass bosons and so our supermassive 
objects should be genuine boson stars. If the phase 
transitions are very recent, the ultra-light bosons 
could belong to the Goldstone sector giving rise to 
mini-boson stars. 

We shall not discuss here any further the for- 
mation processes of boson stars. The reader is re- 
ferred to the recent paper by Mielke and Schunck 
(2000) and references therein. It is apparent that 
for every possible boson mass in the particle spec- 
trum there is a boson star model able to fit the 
galactic center constraints, at least in order of 
magnitude. 

6. Accretion and luminosity 

6.1. Relativistic rotational velocities 

For the static spherically symmetric metric con- 
sidered here circular orbit geodesies obey, 



< = 



,e»-l 8tt 



= e 



M(r) 8tt 



m 



-p r r 
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(47) 
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These curves increase up to a maximum followed 
by a Keplerian decrease. Liddle and Schunck 
(1997) found that the possible rotation veloci- 
ties circulating within the gravitational boson star 
potential are quite remarkable: their maximum 
reaches more than one-third of the velocity of light. 
Schunck and Torres (2000) proved that these high 
velocities are quite independent of the particu- 
lar form of the self-interaction and are usually 
found in general models of boson stars. For in- 
stance, for A = 0, 300 of the Colpi et al. 's 
(1986) choice, U cosh = am 2 [cosh([3 ^/\ : ^p\ 2 ~) - l] , 
and U cxp — am 2 [exp(/3 2 |-0| 2 ) — l] , the maximal ve- 
locities are: 122990 km/s at x = 20.1 for A = 300, 
102073 km/s at x = 4.1 for A = 0, 104685 km/s 
at x = 4.2 for U exp , and 102459 km/s at x = 5.9 
for [/cosh (Schunck and Torres 2000). 

With such high velocities, the matter possesses 
an impressive kinetic energy, of about 6% of the 
rest mass; i.e. to obtain the required luminosity we 
would need that about 10~ 8 to 10~ 7 solar masses 
per year be transformed into radiation. Note that 
the required matter-radiation transfer is at least 
two orders of magnitude smaller than the accretion 
rate towards Sgr A* . 

The maximum rotational speed is attained well 
outside the physical radius of the star, as can be 
seen by computing the dimensionless x value for 
the star radius (e.g. it happens between x ~ 5 
and 15 for A going from to 300). It is interesting 
to note also that the dependence of the maximum 
velocity on A is not very critical, and the same 
process can be operative with mini-boson stars. 
The rotational velocity is dependent on the cen- 
tral density, increasing with a higher value of o~q. 
To obtain large rotational velocities, it is needed 
that the central density of the star be highly rela- 
tivistic, for Newtonian solutions velocities are low 
and quite constant over a larger interval. This is 
consistent with the density constraint of the dark 
object in Sgr A*. 

6.2. The black hole danger 

How can one justify that the accretion onto the 
central object -a neutrino ball or a boson star- will 
not create a black hole in its center anyway. Inter- 
stellar gas and stars, while spiraling down towards 
the center of the object, will begin to collide with 
each other, and may glue together at the center, 



what could be the seed for a very massive black 
hole. Even a black hole of small mass can spiral 
inwards, and if it remains at the center, that black 
hole itself could be the seed. On the other hand, 
stellar formation of massive stars would yield, af- 
ter evolution, to a black hole. Then, we need to 
consider whether there is a mechanism that pre- 
vents the formation of a very massive baryonic 
object -leading inevitably to a black hole- in the 
center of the galaxy. Below we provide an answer 
to this question. 

The key aspect to consider here is disruption. A 
star interacting with a massive object can not be 
treated as a point mass when it is close enough to 
the object such that it becomes vulnerable to tidal 
forces. Such effects become important when the 
pericentre r m j n is comparable to the tidal radius 
(Rees 1988), 

,,= 5xl 0"M./*(£) (M^"> cra . (48) 

Here, -R*,© stands for the radius of the star and the 
sun respectively, while M* t o for the masses. M§ is 
the mass of the central object in millions of solar 
masses. r t is the distance from the center object 
at which M/r 3 equals the mean internal energy 
of the passing star. Alternatively, we may com- 
pute the Roche criterion, a comparison between 
the smoothed out mass of the disruptor and the in- 
ternal mean energy of the disruptee (Krolik 1999). 
This catastrophe radius is 

where R c is given in units of the event horizon, 
R g , of the black hole, and the falling star has den- 
sity p*. Only for black holes with masses smaller 
than 1O 8 M0 and for certain low density stars -red 
giants- we can expect that the disruption happens 
outside the event horizon. This is the reason why 
supplying the material at lower densities, with the 
black hole gravity dominating the situation, can 
generate more power. Rees has also given an es- 
timate of how frequently a star enters this zone. 
When star velocities are isotropic, the frequency 
with which a solar-like star pass within a distance 

^*min IS 

-°- < < 3 (io^)(i55i^)( r r)--' 

(50) 
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where iV* is the star density and a the velocity 
distribution. Disruptions are rare events that hap- 
pens once in about 10000 years. 

6.2.1. Scalar stars 

Because of the similar metric potentials, far 
from the center of the non-baryonic star, the ac- 
cretion mechanism will be the same than that op- 
erative in the Schwarzschild case. If a boson star 
is in the center of the galaxy, the characteristics 
of the tidal radius and the timcscalc of disruption 
occurrence will be similar to those of a black hole 
of equal mass. Stars falling inwards will all be 
disrupted after they approach a minimum radius. 
Contrary to black holes of big masses, which can 
swallows stars as a whole and disrupt them be- 
hind their event horizons, boson stars will disrupt 
all stars -most stars outside and a few inside the 
Schwarzschild radius of a black hole of equal mass- 
at everyone sight. 

In the case of black holes, the most recent sim- 
ulations (Ayal et al. 2000) show that up to 75% 
of the mass that once formed the disrupted star 
become unbound. For boson stars, we argue that 
once the star is disrupted to test particles (with 
masses absolutely negligible to that of the central 
object), and because there is no capture orbits, 
all particles follow unbound trajectories. In this 
sense, all material is diverted from the center and 
the formation of a black hole is avoided. 

It could be worth of interest to perform numer- 
ical simulations changing the central object from 
a black hole to a boson star, to see the aftermath 
and the fate of the debris in an actual boson-star- 
generated disruption. This, however, could well 
not be an easy task: In black hole simulations, the 
minimum radius is maintained still far from the 
center (~ 10 Schwarzschild radius), where Newto- 
nian or Post-Newtonian approximations are valid. 
To actually see the difference between a black hole 
and a boson star it could be necessary to attain 
inner values of radii, where the behavior is com- 
pletely relativistic. 

One case could merit further attention: the pos- 
sible spiraling of a black hole of stellar size. This 
case may complicate the situation since a black 
hole can not be disrupted. However, differences in 
masses are so large that it will behave as a test 
particle for the boson central potential, and will 



be also diverted from the center. Moreover, be- 
ing of stellar size, it will be appreciably influenced 
by other intruding stars, also making it to left an 
static position at the boson star center. 

Finally, it is worth of interest to study other 
possible observational consequences of boson star 
disruption. The ejecta now is 100% of the star 
mass, and may possible produce some long term 
effects in the surrounding medium. In black hole 
scenarios, the bound debris will create a flare as 
it accretes. Does the same happen here? While 
we defer this kind of analysis to another work, 
we mention that if boson stars exists, disruption, 
what was once thought of as an inevitable con- 
comitant of black holes, could now happen in non- 
baryonic environments. 

6.2.2. Neutrino balls 

If the center of the galaxy is a neutrino ball, 
one also has to obtain a mechanism that prevents 
the formation of a very massive baryonic object. 
However, we have to expect very crucial differ- 
ences with a boson star case, caused by the fact 
that a neutrino ball is an extended object and that 
the gravitational potential is shallower. The first 
thing to note is that r t is well within the neu- 
trino ball, and then, stars will traverse the exte- 
rior parts of the ball without being disrupted. In 
doing so, however, the central mass that they see 
at the center will be less than the total mass of the 
ball, and at a distance r — r t the mass enclosed is 
negligible. Disruption can not proceed and other 
mechanism have to devised. We note then that 
the observation of disruption processes in the cen- 
ter of a galaxy is then indicative that a neutrino 
ball is not there. 

When the mass enclosed by the neutrino ball 
is small enough (say, O(10 3 )M Q ), the accretion 
disk will be unstable. This happens about 0.1 - 
1 light years from the center. There, stars which 
could actually form at a rate of 1 per hundred 
thousand years (the actual number will depend on 
the mass of the star) will be probably kicked off by 
intruding stars (Viollier 2000). The absence of the 
disruption mechanism makes this problem worse, 
since given an enough amount of time, it is hard 
to think in compelling reasons by which gas and 
stars are expelled from the center. 

In a published paper, Tsiklauri and Viollier 
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(1999) suggested that matter arriving at the cen- 
ter would be diverted in the form of non-radiating 
jets generated by pressure of the inner accretion 
disk. However, it is not clear that gravity attrac- 
tive forces of the spiraling objects will be smaller 
than gas pressure exerted by the disk. Moreover, 
even when the baryonic mass acquired by the cen- 
tral object during the entire age of the universe is 
at least two orders of magnitude smaller than the 
mass of the central neutrino ball, one should be 
worried if this yields to a black hole of that mass, 
since it may be the seed for a further collapse, or 
appreciably influence the dynamics. The mecha- 
nism considered in the previous paragraph seems 
more adequate than this for solving this question, 
if that is finally possible. To us, it is yet an unclear 
issue in the neutrino ball scenario. 

7. Discussion: how to differentiate among 
these models? 

One of the easiest things one may think of is to 
follow the trajectory of a particular star. This has 
been done by Munyaneza et al. (Munyaneza et al. 
2000) in the case of a neutrino ball. The trajec- 
tory of SI, a fast moving star near Sgr A*, offers 
the possibility of distinguishing between a black 
hole and a neutrino condensate, since Newtonian 
orbits deviate from each other by several degrees 
in a period of some years. However, as soon as the 
central object is not so extended, as in the boson 
star case, this technique is useless (in every case 
in which the pericentre is far than the tidal ra- 
dius), and other forms of detecting their possible 
presence have to be devised. 

It has been already noted that X-ray astronomy 
can probe regions very close to the Schwarzschild 
radius. It is only observations in the X-ray band 
that can study the inner accretion disk as close 
to the center as an event horizon would be. Re- 
cent results from the Japanese-US ASCA mission 
have revealed a broadened iron line feature that 
comes from so close to the event horizon that a 
gravitationally redshift is observed. This is 10000 
times closer into the black hole than what can 
be pictured by HST. In particular, Iwasawa et 
al. (1996) claimed that ASCA observations to 
Seyfert 1 galaxy MCG-6-30-15 got data from 1.5 
gravitational radii, and conclude that the peculiar 
line profile suggests that the line-emitting region is 



very close to a central spinning (Kerr) black hole 
where enormous gravitational effects operate. By 
the way, this is stating that a neutrino ball can 
not be the center of that galaxy. However, as was 
already noted (Schunck and Liddle 1997), a boson 
star could well be a possible alternative, and X- 
ray could be used to map out in detail the form 
of the potential well. The NASA Constcllation- 
X (Constellation-X web page 2000) mission, to be 
launched in 2008, is optimized to study the iron K 
line feature discovered by ASCA and, if they are 
there, will determine the black hole mass and spin 
for a large number of systems. Still, Constellation- 
X will provide an indirect measure of the proper- 
ties of the region within a few event horizon radii. 
A definite answer in this sense will probably be 
given by NASA-planned MAXIM mission (Maxim 
web page 2000), a ^i-arcsec X-ray imaging mission, 
that would be able to take direct X-ray pictures 
of regions of the size of a black hole event horizon. 
Both of these space mission will have the ability 
to give us proofs of black hole existence, or to pro- 
vide evidence for more strange objects, like boson 
stars. 

Very recently, Falcke et al. (2000) have noted 
that gravitational lensing observations of very 
large baseline interferometry (VLBI) could give 
the signature to discriminate among these mod- 
els. Falcke et al. assumed that the overall specific 
intensity observed at infinity is an integration of 
the emissivity (taken as independent of the fre- 
quency or falling as r~ 2 ) times the path length 
along geodesies. Defining the apparent bound- 
ary of a black hole as the curve on the sky plane 
which divides a region where geodesies intersects 
the horizon from a region whose geodesies miss 
the horizon, they noted that photons on geodesies 
located within the apparent boundary that can 
still escape to the observer will experience strong 
gravitational redshift and a shorter total path 
length, leading to a smaller integrated emissivity. 
On the contrary, photons just outside the appar- 
ent boundary could orbit the black hole near the 
circular photon radius several times, adding to 
the observed intensity. This is what produces a 
marked deficit of the observed intensity inside the 
apparent boundary, which they refer to as the 
"shadow" of the black hole. The apparent bound- 
ary of the black hole is a circle of radius 27R g 
in the Schwarzschild case, which is much larger 
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than the event horizon due to strong bending of 
light by the black hole. This size is enough to 
consider the imaging of it as a feasible experiment 
for the next generation of mm and sub-mm VLBI. 
While the observation of this shadow would con- 
firm the presence of a single relativistic object, a 
non-detection would be a major problem for the 
current black hole paradigm. 

Concerning the ideas put forward in this work, 
Falcke et al.'s shadow concept is appealing. In 
the case of a boson star, we might expect some 
diminishing of the intensity right in the center, 
this would be provided by the effect on relativistic 
orbits, however, this will not be as pronounced 
as if a black hole is present: for that case, many 
photons are really gone through the horizon and 
this deficit also shows up in the middle. If a boson 
star is there, some photons will traverse it radially, 
and the center region will not be as dark as in the 
black hole case. A careful analysis of Falcke et al.'s 
shadow behavior replacing the central black hole 
with a boson star model would be necessary to get 
any further detail, and eventually an observable 
prediction. 

We also mention that the project ARISE (Ad- 
vanced Radio Intcrfcromctry between Space and 
Earth) is going to use the technique of Space VLBI 
to increase our understanding of black holes and 
their environments. The mission, to be launched 
in 2008, will be based on a 25-meter inflatable 
space radio telescope working between 8 and 86 
GHz (Ulvestad 1999). It will study gravitational 
lenses at resolutions of tens of ^tarcsecs, yielding 
information on the possible existence (and signa- 
tures) of compact objects with masses between 
10 3 M Q and 10 6 M Q . 

Another possible technique for detecting bo- 
son stars from other relativistic objects could be 
gravitational wave measurements (Ryan 1997). If 
a particle with stellar mass is observed to spi- 
ral into a spinning object with a much larger 
mass and a radius comparable to its Schwarzschild 
length, from the emitted gravitational waves, one 
could in principle obtain the lowest multipolc mo- 
ments. The black hole no-hair theorem says that 
all moments are determined by its lowest two, 
the mass and angular momentum (assuming the 
charge equal to zero). Should this not be so, 
the central object would not be a black hole, 
and as far as we know, the only remaining vi- 



able candidate would be a boson star. In ten 
years time, perhaps, a combination of gravita- 
tional wave measurements, better determination 
of stellar motions, and mm and sub-mm VLBI 
techniques could give us a definite picture of the 
single object at the center of our Milky Way. 

From a theoretical point of view, developments 
on gravitational lensing theory in very strong field 
regimes will be of extreme importance since, for 
objects like Sgr A*, the standard weak field theory 
does not hold and new effects have to be expected 
(Virbhadra and Ellis 2000). 

8. Conclusions 

We have shown that boson stars provide the ba- 
sic necessary ingredients to fit dynamical data and 
observed luminosity of the center of the Galaxy. 
They constitute viable alternative candidates for 
the central supermassive object, producing a the- 
oretical curve for the projected stellar velocity 
dispersion consistent with Keplcrian motion, rel- 
ativistic rotational velocities, and having an ex- 
tremely small size. 

Other singularity-free models were as well con- 
sidered, as non-baryonic fermion stars (e.g. neu- 
trino, or gravitino condensations). In this case, 
the object is sustained by its Fermi energy, while 
in the boson star case, it is the Heisenberg's Un- 
certainty Principle which prevents the system from 
collapsing to a singularity. Due to this fact, bo- 
son stars are genuine relativistic objects where a 
strong gravitational field regime holds. 

This difference in the relativistic status of both 
objects is not trivial. While fermion neutrino balls 
are extended objects, boson stars mimic a black 
hole. Disruption processes can not happen in a 
fermion condensation, and it has the unpleasant 
consequence of not providing with a straightfor- 
ward mechanism by which stars could be diverted 
from the center, and through which finally avoid 
the formation of a massive black hole inside the 
condensate. 

The formation of boson stars, neutrino balls, 
and black holes can all be competitive processes. 
Then, it might well be that even if we discover that 
a black hole is in the center of the galaxy, others 
galaxies could harbor non-baryonic centers. In the 
case of boson stars, only after the discovery of the 
boson mass spectrum we shall be in position to de- 
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termine a priori which galaxies could be modeled 
by such a center. Observations of galactic centers 
could then suggest the existence of boson scalars 
much before than their discovery in particle physi- 
cists labs. 
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A. Appendix 

For reader's convenience, we quote here the dimensional conversion for the radius and the mass of a boson 
star. Using the value of 1 GeV in cm -1 , and taking into account the dimensionless parameter x = mr, we 

r ^ = M&vj 6 - 38 x 10 " 33 - < A1 ) 

For the mass, recalling that M = M(x)mp l /m, we get 

M[1O 6 M ] = ^ X } T , 1.33 x 1(T 25 . (A2) 
m[GeVJ 

In the case where A>1, both right hand sides of the previous formulae get multiplied by A 1 / 2 . 
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Fig. 1 . — Enclosed mass in the center of the galaxy 
together with observational data points. See dis- 
cussion in the main text. 

Fig. 2. — Projected velocity dispersions: observa- 
tional data and fit using a boson star model. See 
discussion in the main text. 

Fig. 3. — Mass and number of particles, in dimen- 
sionless units, for a A = 10 boson star (Colpi et 
al. (1986)) potential. 

Fig. 4. — Stability analysis for A = 10 configura- 
tions. Catastrophe theory ensures us that the first 
branch, which includes the (0,0) point, is the only 
stable one. 

Fig. 5. — Boson star mass as a function of x, for a 
A = 0, (t(0) = 0.1 model. We show a Boltzmann- 
like fitting (dotted curve, but almost everywhere 
on the solid curve) and its residue (solid lower 
line). 

Fig. 6. — Boson star metric potentials g rr = 
(dashed) and gtt — e v (solid). Boson star param- 
eters are A = 0, cr(0) = 0.1. 

Fig. 7. — Boson (solid) and Schwarzschild 
(dashed) effective potentials for massless particles, 
b = I 2 /E 2 ^. The mass for both, the black hole and 
the boson star, was taken as M — 0.62089mp[/m. 
The maximum in the black hole case happens, in- 
dependently of I, for r/M = 3. 

Fig. 8. — Boson (solid) and Schwarzschild 
(dashed) effective potentials in the case of mas- 
sive particles. The mass of the central object is 
as in the previous figure. The three curves corre- 
spond to l 2 /M 2 = 0, 12 and 15 (from bottom to 
top). 

Fig. 9. — Constraint in the boson star fundamen- 
tal parameters which gives rise to an object of two 
million solar masses within approximately ten so- 
lar radius, consistent with the mass of the central 
object in our galaxy. 

Fig. 10. — Constraint in the non-topological soli- 
ton fundamental parameters which gives rise to an 
object consistent with the mass of the central ob- 
ject in our galaxy. It is specially marked the usual 
case in which $ = m. 
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